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Figure captions 

Figure 1 The momentum dependence of the ratios R(r,p,N) (3.1) at r = |,|,|,| 

and N = 4, 6, 28. The solid curve is the prediction of reference [12]. 
Figure 2 The finite size behaviour (3.9) of p(t,p,N) - defined in (3.2): (a) The N~ 2 

dependence at p/ir = |,|,| and t = 1,5,10. (b) The coefficient c(t,p) versus t 

(3.3) at fixed p/ir = \, |, ±, §, f . 
Figure 3 The estimate p(t,p) of the thermodynamical limit versus t at fixed momenta 

p/ir = j, |, |, |, |, |, |, |, |, |, |. The inset shows the magnification for < t < 1. 
Figure 4 Comparison of the estimate p(t,p) with: (a) the ansatz of Miiller et al. [12] 

(b) the modified ansatz (4.3). 
Figure 5 The momentum dependence of the parameters in the modified ansatz (4.3): 

(a) the modification u(p) of the high-frequency cut-off (4.4) (b) The amplitude 

A(p). 
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K^ n \p) listed in section 4. Deviations from our estimate for the thermodynamical 
limit appear for larger momenta and euclidean times r. 
• In the second ansatz (4.3) the momentum dependences of the amplitude A(p) and 
of the modification u(p) for the high-frequency cut-off (4.4) have been chosen such 
that the spectral moments K^ n \p) with n = —1, 0, 1 are reproduced. In particular, 
the zero momentum behaviour of K^~ l \p) and K^ l \p) enforces the vanishing of 
the amplitude A(p = 0) and a well defined modification (4.5) of the high-frequency 
cut-off (4.4). The second ansatz is in better agreement with our estimate for the 
thermodynamical limit of the dynamical structure factor. 

We are aware of the fact, that the relevant excitation spectrum in the spin-1/2 
antiferromagnetic Heisenberg model is unbounded. One might ask how to interpret 
the 'effective' high-frequency cut-off (4.4) and the fact that it is above the two-spinon 
cut-off for p < 7r/3 and below for p > 7r/3? Adding in the dynamical structure factors 
high-frequency excitations (u > Ll> 2 (p)) - which are not of the two-spinon type - will 
always lead to an increase of the effective high-frequency cut-off (4.4). Therefore, such 
contributions are responsible for the increase of the effective high-frequency cut-off (4.4) 
for p < 7r/3. On the other hand they cannot explain the lowering of the effective high- 
frequency cut-off for p > 7r/3. 
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This yields a prediction for the static structure factor in the low-momentum limit: 

K W( P ) ^ 1.084... • (4.6) 

47T 

which is in good agreement with the results of the finite size analysis [17]. For 
7r/14 < p < 137T / 14 , we fix u(p) in such a way that the ratio (3.6) of the moments with 
n = 1 and n = is reproduced correctly. u(p) is shown in figure 5(a). It turns out to be 
smaller than 1 for p > 7r/3 and larger than 1 for p < 7r/3. The extrapolation to p = 
seems to meet the value given in (4.5). Finally, we determine the p-dependence of the 
amplitude A(p) from the moment sum rule (3.8) with n = 1. A(p) is shown in figure 
5(b). The extrapolation to p = is compatible with the behaviour given in (4.5). The 
values of A(p) and u(p) at p = ir: 

A(tt) = 1.88(8), u(tt) = 0.63(5), (4.7) 

have been obtained in reference [17] from a fit to the static structure factor of the form: 

K {*\ p ) ^-^l\ n (i-P) + const . (4.8) 

In figure 4(b) we compare the estimate for p(t,p) with the prediction pMM(t } p) } 
which follows from the ansatz (4.3) with the momentum dependent amplitude A(p) and 
the modified high-frequency cut-off (4.4). The difference: 

A MM (t,p) = p MM (t,p) - p(t,p) (4.9) 

is obviously much smaller than (4.1), which is shown in figure 4(a). On the other hand 
(4.9) is not yet zero. In particular for large t and p = 37r/4 we still observe sizeable 
deviations form our estimate p(t,p) for the thermodynamical limit, leaving room for 
further improvements of the ansatz (4.3). 

5. Conclusions 

The euclidean time representation (1.8) is particularly suited for a study of finite 
size effects in dynamical structure factors. In this paper we investigated the spin- 
1/2 antiferromagnetic Heisenberg model (1.1). We found that finite size effects in the 
euclidean time representation die out with N~ 2 for those momentum values p where 
there is a gap (1.4) between the groundstate and the first excited state with momentum 
p. We then compared our estimate for the thermodynamical limit with two predictions. 
Each of them follows from an ansatz of the type (4.3) for the spectral (u)- representation 
of the dynamical structure factor: 

• The first ansatz - i.e. (4.3) with momentum independent amplitude A(p) = A and 
u(p) = 1 - has been proposed by Miiller et al. [12]. It takes into account two-spinon 
contributions only and respects the properties (i) and (ii) on the frequency moments 
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The error in our estimate of p(t,p) is presumably less than the width of the lines in 
figure 3. This optimistic view appears to be justified by the clean signal for the finite 
size behaviour (3.9) seen in the data. 

4. A modified ansatz for the dynamical structure factors 

In figure 4(a) we compare our estimate for p(t,p) with the prediction pM(t } p) } which 
follows from the ansatz (1.6). The difference: 

A M (t,p) = PM (t,p)-p(t,p) (4.1) 

is negative for p < 7r/3 and positive for p > 7r/3. It increases with t and p for p > 7r/3. 
The increase of the deviations with p has been observed also in the static structure factors 
[17]. Before presenting modifications of the ansatz, which are in better agreement with 
the estimate of the thermodynamical limit let us briefly review the properties of the 
frequency moments K^ n \p) - defined in (3.7): 

(i) The moment (3.8) with n = 1 is determined by the groundstate energy. 

(ii) The higher moments with n odd are known [13] to be polynomials of order n in 
cos p. 

(iii) The moment with n = — 1 is related to the static susceptibility x{p)i which yields 
in the zero momentum limit [12]: 

K^\p) ^1* = i (4-2) 

(iv) The moment with n = is identical with the static structure factor, which has been 
determined in reference [17] from a finite size analysis in the momentum interval 
tt/14 < p < 13tt/14. 

The ansatz of Miiller et al. [12] respects the first two properties. However, the first and 
third one cannot be satisfied simultaneously with this ansatz. We therefore would like 
to propose a modification here: 

s(u,p) = n ^£!=e(u-u 1 )e{Lu 2 -u) (4.3) 

with a momentum dependent amplitude A(p) and a modified high-frequency cut-off: 

u 2 {p) = u(p)u 2 (p). (4.4) 

Properties (i) and (iii) are satisfied now if we choose: 

u{p) 1.274..., A(p) ^ 0.7475... • p. (4.5) 



follows from a Taylor expansion of (1.8) around r = 0. The coefficient: 



T= o K(0) (P) 

is given by the first two moments of the dynamical structure factor: 



K {n \p) = / duu n S(u;,p). 



(3.6) 



(3.7) 



K^°\p) is just the static structure factor. It has been determined numerically for 
7r/14 < p < 137T / 14 in reference [17]. The first moment is known analytically [13]: 

K {1) (p) = (1 -cosp)e (3.8) 

where e = (4 In 2 — l)/3 is related to the groundstate energy. The coefficient P2(p) is 
given in table 3. 



Table 3. The ratio (3.6) of frequency moments K^/K^ versus p. 



p/tt 


1/14 


1/13 


1/12 


1/11 


1/10 


1/9 


1/8 


1/7 


P2(P) 


0.7373 


0.7919 


0.8552 


0.9296 


1.0180 


1.1249 


1.2566 


1.4226 


P/tt 


1/6 


1/5 


1/4 


2/7 


1/3 


2/5 


3/7 


1/2 


P2(P) 


1.6368 


1.9228 


2.3173 


2.5705 


2.8662 


3.1930 


3.2997 


3.4741 


P/tt 


4/7 


3/5 


2/3 


5/7 


3/4 


4/5 


5/6 




P2(P) 


3.5133 


3.4912 


3.3579 


3.1933 


3.0353 


2.7634 


2.5496 




p/tt 


6/7 


7/8 


8/9 


9/10 


10/11 


11/12 


12/13 


13/14 


P2(P) 


2.3805 


2.2432 


2.1301 


2.0353 


1.9545 


1.8847 


1.8237 


1.7698 



For all values of p and N we observe a very smooth and almost linear behaviour 
of p in t, except near t = 0, where the quadratic behaviour (3.5) is visible. Finite size 
effects can be described with high accuracy by: 

p(t,p,N) = p(t,p,N = 00) + (3.9) 

as is demonstrated in figure 2(a) for p/ir = |,|,| and t = 1,5,10. The coefficient 
c(t,p) increases with t and with p for p > 7r/3, as can be seen in figure 2(b). The 
resulting estimates p(t,p) = p(t } p } N = 00) for the thermodynamical limit are plotted 
in figure 3 for p/ir = j , |, |, |, |, |, |, |, |, |, |. Here we have included momentum values 
p = kir /7, k = 1, 2, 6, which are realized on two systems (N = 14, 28) only - assuming 
that the finite size behaviour is of the form (3.9). Note that p(t,p) is strictly monotonic 
in p and t. In order to reach this property we have changed the definition of the variable 
t. (3.3) differs from the definiton given in reference [15] by the factor Jui(p). 
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Table 2. Comparison of leading excitation energies and transition probabilities with 
numerically exact results. 

N u>i(tt) wi(tt) 

18 0.4824988997812 0.785322982816 

20 0.43589108737 0.771429257660 

22 0.3975468223202 0.7592200871420 

24 0.365442071611 0.748370457307 

26 0.3381648484657 0.73863676355 

28 0.314699900401 0.729832566394 



dependence of the ratios: 

b(0,p,N) 

for t = |,|,|,| fixed and N = 4, 6..., 28. The data points for N > 8 scale in the 
momentum p except for p = ir. Therefore the resulting curves at fixed r represent 
already the thermodynamical limit R(r } p } N = oo) for p < ir. The curves start with 
R = 1 at p = 0, have a minimum at p = Po(t) ir/2 and approach again R = 1 for 
p = it. The latter is a consequence of the nonintegrable infrared singularity in (1.5). The 
solid curves in figure 1 represent the prediction of the ansatz (1.6) of Miiller et al. [15]. 
This prediction is in good agreement with our finite system results for small p values. 
Deviations emerge for larger p values, which will be analyzed in the next section. 

For increasing values of r finite size effects increase. They become visible in the 
quantities p(t } p } N) } which are related to the ratios (3.1) via: 

R,T -"- N) = TTW^N)- (3 ' 2) 

In terms of the variable: 

t = Y^(p)r exp[cui(p)r] (3.3) 
the large t behaviour: 

p(t,p,N = oo)^t (3.4) 



is linear. (3.4) originates from the threshold singularity (1.6) which is projected out in 
(1.8) in the large r limit. 
The small t behaviour: 

p(t,p,N = oo)^^ 2 (3.5) 
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Table 1. Comparison of leading excitation energies and transition probabilities with 
the results of complete diagonalization [15] on a ring with 16 sites. 



S(r = 0, p = tt/4) = 2.982766323178 -10" 1 S(r = 0, p = tt/2) = 6.794375761266 -10" 



w n (7r/4) 


w n (7r/4) 






w n (ir/2) 








o.oouuuioooooyo 


8.603551483282 


•10 


1 


4.29305867857 


1.0940110661 


•10" 3 


4.197135363571 


1.288323108707 


•10" 


1 


4.55019377043 


8.579155788 


•10" 3 


4.59757074255 


6.12984208 


•10" 


3 


5.00373864011 


7.2716867 


•10" 4 


4.73871528573 


1.48596359 


•10" 


3 


5.568925 


2.792 


•10" 4 


4.902809151 


9.727121 


•10" 


5 


5.781 


4.38 


•10" 4 


5.35132697 


1.723352 


•10" 


4 


5.8450 


1.66 


•10" 3 


5.8903 


3.01 


•10" 


4 


6. 




O(10" 4 ) 


5.9 


1.875 


•10" 


3 


6. 




O(10" 4 ) 


6.47 


8. 


•10" 


5 



S(r = 0, p = 3tt/4) = 1.3230534343023 S(r = 0, p = tt) = 4.292303508279 



w„(3tt/4) 


w„(37r/4) 








w„(ir) 


2.6381304345681 


7.717163620986 


•10" 


1 


0.540379364500 


8.01345217378 


•10" 1 


3.411532042825 


7.7506497226 


•10" 


4 


2.79206117219 


1.4109158151475 


•10" 1 


4.330655742114 


1.8574615465811 


•10" 


1 


4.668596605332 


4.6691968837931 


•lO" 2 


5.03583236201 


8.985257272 


•10" 


4 


5.475947091450 


1.636197084 


•10" 4 


5.440678 


4.4375 


•10" 


4 


5.9070165813 


1.044901942 


•lO" 2 


5.46270297 


3.950172 


•10" 


2 


5.994081 


1.03778 


•10" 6 


5.980403 


8.0786 


•10" 


5 


6.573253 


1.55186 


•10" 4 


6.44 


2.3 


•10" 


5 


6.80283 


6.702 


•10" 5 


6.643 


5. 


•10" 


4 


7.1 




0(10" 



the recursion method with L = 40 iterations is given in table 2. Again we only list 
those digits which agree with the numerically exact result. The latter is obtained from 
a determination of the groundstate in the channels with total spin S = and S = 1, 
respectively. 



3. Finite size analysis of the dynamical structure factors in the euclidean 
time representation 

For noncritical momenta p < tt and euclidean times r not too large finite size effects 
are small in the dynamical structure factors (1.8). In figure 1 we show the momentum 
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with matrix elements: 



</*'|0|/*> 



' 1 : 


k' = k + 1 


a k : 


k' = k 




k' = k-l 


, : 


else 



(2.6) 



in the basis \fk}- The eigenvalues of yield the excitation energies u = E n — E and 
the eigenvectors represent the excited states \n) of (H — E ) in the basis \fk}'- 

W = t\h)j7Try (2-7) 

The zero components (fo\n) of the eigenvectors and the excitation energies E n — E 
determine the dynamical structure factor: 

S(T,p,N) = J2\(Mn)\ 2 exp[-T(E n -E )]. (2.8) 

n 

Performing the iteration (2.2) numerically, one finds that the orthogonality of the vectors 
\fk) is lost after a certain number of steps due to rounding errors. On the other hand 
it turns out that the energies E n — E and transition probabilities of the low-lying 
excitations with n < 10 can be obtained already with high accuracy by truncating the 
problem. We stop the iteration after L steps and diagonalize the truncated Lx L matrix 
O. As an example we show in table 1 the first 9 excitation energies and normalized 
transition probabilities 

Mp , N) jM^mi ( , 9) 

S{0,P, N) 

on a ring with N = 16 sites and momenta p/ir = |, |, |, 1. 

The approximate solution given in table 1 is obtaind by truncating the iteration 
after L = 40 steps. We only list in table 1 those digits which agree with the numerically 
exact solution given in [15]. In other words, the number of quoted digits measures the 
accuracy of the approximation. Note that the energies and transition probabilities of the 
first four excitations are reproduced correctly with 8 and more digits. The accuracy for 
the remaining 5 excitations is less impressive. However, this inaccuracy has practically 
no effect on our evaluation of the dynamical structure factors in the euclidean time 
representation (1.8), since the contributions of the higher excitations are suppressed 
twofold. In addition to the exponential damping factor exp(— lot) in (1.8) it turns 
out, that the transition probabilities themselves drop rapidly with u. Therefore we 
expect, that the determination of the excitation energies and transition probabilities 
by means of the recursion method will yield as well reliable results for larger systems 
with N = 18, 20, 28. In these systems we have checked the energy and the transition 
probability of the first excitation with momentum p = ir. The result obtained with 



3 



in u. In order to extract the thermodynamical limit it is useful [15] to consider the 
Laplace transform of (1.2): 

oo 

S(t, P ,N)= J du S(u,p,N)exp(-uT) (1.8) 

= (0\st(p)exp[-r(H - E o )]s 3 (p)\0). 

It can be interpreted as a euclidean time (r)-representation of the dynamical structure 
factor. 

In this paper we want to propose a new method to compute the dynamical structure 
factors. It is based on the recursion used in [16] as input for the continued fraction 
approach. In section 2 we will demonstrate, how the energies and transition probabilities 
for the low excitations can be obtained directly from the recursion approach. In section 
3 we will present our results on the dynamical structure factors (1.8) for systems 
with N = 4, 6, 8, 28. Finite size effects will be analyzed and an estimate for the 
thermodynamical limit will be given. This estimate is compared in section 4 with the 
prediction (1.6) of Miiller et ai. We also propose a modification of this ansatz, which 
yields better agreement with the estimate of the dynamical structure factors in the 
thermodynamical limit. 



2. The Recursion method 



Following [16], (1.8) can be computed by iteration. For this purpose one expands the 
euclidean time evolution of the 'initial' state |/ ) = s 3 (p)\0): 

oo 

exp[-r(tf - E o )]s 3 (p)\0) = £ D k (r)\f k ) (2.1) 

k=0 

in terms of an orthogonal basis \f k ) which is constructed recursively by application of 
the Hamilton operator: 

(H - E )\f k ) = \f k+1 ) + a k \f k ) + k = 0, 1, 2, ... (2.2) 

(fk\H - E \f k ) 

(fk\fk) ' /s = ' 1 ' 2 '- t 2 - 3 ) 

,2 _ (fk\fk) J. _ 1 9 (O A\ 

6 *-</*-ll/*-l> '- 1 ' 2 '-' (2 - 4) 

with = and 6q = 0. On finite systems the iteration (2.2) will terminate 

after L steps, where L is the dimension of the Hilbert space spanned by the states 
(H — Eq) 1 s 3 (p) |0), / = 0,1, L — 1. The iteration generates a tridiagonal L X L-matrix 
O: 



(H - E )\f) = 0\f) 



(2.5) 
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1. Introduction 



Quantum spin systems with known dynamical behaviour are rare. The spin 1/2 
XX-model can be mapped on a free fermion system [1] and the dynamical spin-spin 
correlators can be computed analytically [2-5]. A second example is the Haldane- 
Shastry model [6,7] - an isotropic Heisenberg model with couplings which decrease 
with the inverse square of the distance between two spin operators. In this model only 
two-spinon excited states contribute to the zero temperature dynamical structure factor 
[8,9]. In this paper we are concerned with the familiar antiferromagnetic Heisenberg 
model with nearest neighbour coupling: 



N 



H = 2 J2 s(x)s(x + 1) (1.1) 

x=l 

and periodic boundary conditions. The dynamical structure factors at T = 0: 

S(u,p, N) = E% " (En - E o ))\(n\s 3 (p)\0)\ 2 (1.2) 

n 

are defined by the transition probabilities |(n|s 3 (p)|0)| 2 from the groundstate |0) to the 
excited state \n) with energies E and E n , respectively. The transition operator 

1 N 

S3 ^ = ~7m ^ exp(ipx)s 3 (x) (1.3) 

is just the Fouriertransform of the spin-operator s 3 (x) at site x. The lower bound for 
the excitation energies 

Ll>i = iOi(p) = irs'mp (1-4) 

has been computed by Cloizeaux and Pearson [10]. The infrared behaviour at p = ir: 

S(tu^p = 7r, N = oo) t '^ to~ l (1-5) 

has been investigated in reference [11]. G. Miiller and collaborators [12] proposed an 
ansatz, which takes into account only two spinon excitations: 

S{u,p) = A . (1.6) 

The upper bound of these excitations 

lo 2 {p) = 27rsin| (1.7) 

has been obtained from the Bethe ansatz. The ansatz (1.6) has been shown [13] to 
satisfy certain conditions on the spectral moments for the dynamical structure factor. 
Moreover, it was successfully applied to the description of neutron scattering data [14]. 
On finite systems, the dynamical structure factors are defined as (^-function contributions 
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